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ABSTRACT 

We consider the joint probability distribution function for the mass contrast 
and angular momentum of over-density regions on the protogalactic scale and in- 
vestigate the formation of massive black holes at redshift 10. We estimate the 
growth rate of the angular momentum by the linear perturbation theory and the 
decay rate by the Compton drag and apply the Press-Schechter theory to obtain 
the formation rate of massive black holes, assuming the full reionization of the 
universe at z = z. ^> 10. We find the correlation between the mass contrast 

ion 

and angular momentum vanishes in the linear theory. However, application of the 
Press-Schechter theory introduces a correlation between the mass contrast and an- 
gular momentum of bound objects. Using thus obtained probability distribution, 
we calculate the mass fraction of black holes with M ~ 10 6 - 10 8 M Q in the universe. 
We find that it crucially depends on the reionization epoch z. . Specifically, for 
the standard CDM power spectrum with the COBE normalization, the condition 
z. ~ 500 must be satisfied to reproduce the observed number density of QSOs. 



1. Introduction 



Now in the cosmology, one of the least known periods lies between the recom- 
bination time and the redshift z ~ 5 at which most distant QSOs are observed. 
This era has given us very little observational information, but is believed to be 
a stage during which various important events and processes for the structure for- 
mation occurred. Among them, one that has recently attracted much attention is 
the reionization of the universe. If a sufficient number of massive stars had formed 
soon after recombination, the universe would be reionized at very high redshifts 
z £ 100. Although we have no firm theoretical ground to estimate the number 
of such massive stars at the moment, there are various observational reasons which 
favor the reionization scenario. 

Concerning the effect of reionization on the large scale structure, anisotropics 
of the cosmic microwave background (CMB) at intermediate angular scales (~ 
arcminute to degree scales) would be reduced considerably to loosen the constraint 
on the amplitude of the primordial density perturbation on corresponding scales 
(~1 to 100 Mpc). [21 In particular, this may help explaining the existence of large 
scale velocity fields without contradicting the observed upper bound of the degree 
scale CMB anisotropyP 1 

Reionization of the matter also strongly affects the formation of the small scale 
structure of the universe. Right after the recombination, the Jeans mass for baryons 
is ~ 10 6 M Q and in the absence of reionization it gradually decreases as the uni- 
verse expands. However, if the universe was reionized, there would have been a 
considerable increase in the Jeans mass. This would mean the interruption of the 
small scale structure formation. This situation can be avoided if the dominant 
component of the universe is some cold dark matter (CDM), since the Jean mass 
for CDM is negligibly small. Then the baryons fall into the CDM potential well 
to form protogalactic clouds and they may start collapsing once the optical depth 
becomes sufficiently small and the cooling time is short compared with the free-fall 
time. Once we assume such a scenario, then a matter of interest is whether those 
protogalactic clouds would collapse completely to form black holes, since the for- 
mation of such massive black holes at high redshifts can explain the existence of 
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QSOs and/or AGNs at z ~ 5. 

Closely related subjects have been discussed previously by various authors. It 
is already shown that the number density of QSOs can be explained [4,51 in the 
context of the Press-Schechter theory | 6 ' 71 which assumes a simple extrapolation 
of linear theory to the non-linear growth of density perturbations. However, they 

[^] 

neglected the effect of angular momentum. As was pointed out by Peebles, the 
angular momentum gained during the growth of the density perturbation may be 
large enough to support the cloud against runaway collapse into a black hole. ^ 
Loeb' 131 and Umemura et al!" 1 calculated the angular momentum barrier of the 
density perturbation and showed that the angular momentum of the cloud can 
be reduced sufficiently to form a black hole by Compton drag, i.e., the inverse 
Compton scattering with CMB photons, provided that the reionization epoch was 
early enough. Then what should be done now is a statistical analysis to evaluate 
the fraction of baryons which turns into black holes, given a primordial density 
perturbation spectrum. 

In this paper, assuming the reionization scenario as mentioned above, we inves- 
tigate the formation of massive black holes (M ~ I0 6 - 10 8 M Q ) at high redshifts 
and calculate the mass fraction of those black holes in the universe to test whether 
such a reionization scenario can explain the observed number density of QSOs. We 

[15] 

assume a CDM dominated universe with the Harrison- Zeldovich spectrum and 
normalize the amplitude by the COBE-DMR result at 10° J 161 Reionization is as- 
sumed to have occurred at once at z = z. , which we regard as a free parameter. 

ion 1 ° 1 

The paper is organized as follows. In section 2, we estimate the angular momentum 
barrier of a cloud by linear theory and calculate the joint probability distribution 

— * — * 

P( <*>(r ); t) for the angular momentum J{r' ) and mass contrast <5(r ) of pro- 
togalactic clouds. We find there is no correlation between the angular momentum 
and mass contrast in linear theory. In section 3, we apply a Press-Schechter-like 

— * 

theory to the joint probability distribution function P( J(r Q ), <5(r Q ); t) to formulate 

f Although there have been arguments on the proper definition of the angular momentum, 19 121 
in this paper, we adopt the one according to Peebles.' 81 See section 2 below for further 
discussion on this point. 
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a method to calculate the mass fraction of black holes in the universe, taking into 
account the loss of the angular momentum by the Compton drag. In doing so, a 
correlation between the angular momentum and mass contrast is introduced. In 
section 4, we show our results of calculations. Within the context of CDM scenario 
with the COBE normalization, we find z. must be greater than 500 to explain 

' ion ° 1 

the observed number density of QSOs. Finally, section 5 is devoted to summary 
and discussion. 



2. Linear theory 

In this section, we derive the probability distribution function for the mass 

— * 

contrast 8(r ) and the angular momentum J(t~ ) averaged over a sphere of comoving 
radius r Q in the linear perturbation theory. 

For definiteness, we adopt the standard picture of the very early universe, 
namely a cosmological model based on the inflationary scenario, and assume = 1. 
As the epoch we are interested in is the matter dominated era, we set 

I x 2/3 

where p Q is the cosmic mean density at the present time and t Q is the age of the 
universe. Furthermore we assume the universe is cold dark matter dominated and 
baryons contribute only a small fraction of the cosmic density, O fe <C 1. Thus we 
assume £ = 8 at the stage of our interest. 

[^] 

Following Peebles, we define the angular momentum measured not from the 
geometrical center of the sphere but from its center of mass. Hence, 

«('-<,) : = 5^3 / Ht,r)d 3 r, 

V 'f (2.2) 

J(r ):=O fe a 4 / p(t)(l + 6(t, r))(r - R) x u(t, r)d 3 r, 

\r\<r 

where u(t,r) is the peculiar velocity field, r represents the comoving coordinates, 



-A- 



and R is the center of mass coornidate of the sphere, 



R= 7 / 8(t,r)rd 3 r. (2.3) 

47rrj] J 

\r\<r 

Note that our definition of the angular momentum is second order in the perturba- 
tion amplitude, which is different from the first order expression one obtains in the 
Lagrangian picture as in some other recent works! 9 121 But provided we consider 
the angular momentum gain until 8 becomes of order unity, the resulting magni- 
tude of the angular momentum is not so different from each other. Furthermore, in 
our models discussed in section 3, if the angular momentum loss due to Compton 
drag is absent, the resulting spin parameter A := J \fEG~ 1 M -5 / 2 turns out to be 
mass-scale independent. Here E denotes the binding energy of the density pertur- 
bation. Specifically, we obtain A ~ 0.025 for Model 1 and A ~ 0.08 for Model 2 (see 
Eqs.(3.f ) and (3.2) below), while the value A ~ 0.05 has been obtained in numerical 
simulations.' 171 Hence after all, our definition is not a bad approximation for the 
primodial angular momentum. Note also that, as we are interested in the baryonic 
component, the angular momentum we have defined is that of baryons. 

Under the assumption of the rotation-free velocity field and the Gaussian nature 
of the primordial density perturbations, which are also naturally predicted in the 
inflationary universe scenario, the density perturbation 8 and the velocity field 
u(t, r) may be expressed in terms of the Fourier components as 

8(t,r) = j dka h (t)z k e^\ 

2z ci k , ik-v 



/ \ f 7i 2« a k 

u(t,r) = / dk —z(T,{t)z,e 

y ' ; J 3 t k 2 kV J k 



where cr^it) is the square root of the density power spectrum which is proportional 

to the scale factor in the present case, and is an independent stochastic variable 

for each k except that z ^ must be identified with z because of the reality of 8. 

The distribution of z, is normalized as 
k 

/ x dx.dy, I x, 2 -\-y, 2 \ 
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where 



z k = 71 ( x k + i y k ) ' x - k = x k> y-k = -y k - ^ 

In terms of the these Fourier variables, Eq.(2.2) is rewritten as 

^ r o) = Eva 



k 



kk' 



(2.7) 



where 

h :=W M( kr oK> 

8vr„ (ar n ) 5 kxk' 
9 — ; 

( A? ^ A? , A? ^ A? 



(2.8) 

and W ' M and are the Fourier transforms of the top hat window function and 
its derivative, respectively, which are given by 

W M (y):=3^, IP>):=15^. (2.9) 

Hereafter we drop the argument r Q , for notational simplicity. Since the sta- 
tistical nature of the density distribution is all known, we can calculate the joint 
probability distribution of the two variables by performing the following integration, 



P M [J.6\t 



J V kk' / 



k 



xp z{ z k)Ti dx k d yk> 
k 

where the primed product Fj' represents a product over a half of the Fourier space 
k for which and are independent, e.g., those which satisfy k x > 0. Using the 



- 6 - 



Fourier transform of Dirac's ^-function, we rewrite Eq.(2.10) as 



P js[ J > S ' t ) = l I -^3^exp(zn<5)exp (is ■ J 



k \ kk' / 

^Tl x kh ex P \-\ E' M fcfc'ffcffc 

fc / V kk' / 



(2.111 



k \ kk 

exp — in 



where 



(2.m 



M jk* ■= s k k . + 2 *fk k . ■ * , 

M, , = 8,, - 2if- ■ s, 

kk' kk' J kk' ' 
fkk' := fkk' + ^fc-fc' ' 

— * — * — * 

fkk' := fkk' ~ fk-k' ' 

and Y^,' i s defined similarly as Fj'. After the integration with respect to dx^ , dy^ 
and c/n, we obtain 



P J8 (J,8;t) = I exp (z3 • J) (det (mm))" 1 

/ ,2 \ (2-13) 



, ■ exp . , 

4vr £' M~l g k g k , V 4 E kk , M£g k g k , , 



fcfc' kk'^k^k 



We see that from this expression the possible correlation between 8 and J arises 
only through the ^-dependence of the last exponential factor, 



/ 8 

exp 



4 V' Mr}g,g,, 
^ kk' kk' y k y k' 



Interestingly, this term turns out to be ^-independent. Consequently there exists 

— * 

no correlation between 8 and J in the linear perturbation theory. Now let us show 
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this fact, or that 



(2.14) 



is ^-independent. To show this, we expand C in the power series of s as 

oo 

C =Y, C N> ( 2 - 15 ) 



N=0 



where 



k 



c n= E hh>(-™- fkkj-i-™- / fej V -i fe,; " 



fc,fclfc2..fcjV-l,fc' 



(2.16) 



^From the definition of given in Eq.(2.12), we can rewrite the primed sum over 
the last index k' in C N to the unprimed sum as 

c n= E' ^ibOE^c-^-V^')- ( 2 - 17 ) 

fc,fcifc 2 ...fcjv-i fc' 



Now, since / given in Eq.(2.8) has the form, 

fc_ zv _ 1 fc 



^ fe/ = k x fc' (W(|fc + k'\,k, k') -W(\k + k'\,k', k) , (2.18) 

it is easy to see that for each k there exists one and the only one k whose contri- 
bution exactly cancels that from k (see Fig.l). Thus C N = for all N > 1 and we 
obtain 

C = E'^ = ^E^- ( 2 - 19 ) 

As a result, the joint probability distribution factorizes: 

P J6 (J,6;t) = P 6 (6;t)-Pj(J;t), (2.20) 
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where 



PJS;t) 



Pj(J;t) 



d 3 s 



exp 



'2(< 2 ) 



9k 



(2.21) 



(2tt) 



exp lis - J) (det I MM) 



W 2 



In contrast to the apparent Gaussian nature of P , the explicit form of Pj is 
hard to obtain. One case that can be analytically calculable is the limit of large 

— * 

| J |, for which we may approximate the determinant factor by exponetiating it, 



(det MM} 



-1/2 



exp 



iTrln(MM) , amd expand the exponent in powers of s 



to O(s^). The resulting form of the probability distribution is Gaussian, 



PJJ;t)d 3 J 



1 



3/2 



exp 



2(J 2 )/3 



d 6 J\ J 



kk' 

(2.22) 



Although there is no justification for this to hold when \ J\ is small, hereafter, we 
assume the above Gaussian form for Pj} 



3. The collapsed mass fraction formula 

Now we have the approximate factorized joint probability distribution, valid at 
the linear stage. Its temporal behavior is completely determined by that of (8 2 ^ 
and (^J* 2 ^ ■ F° r (3 2 )i we simply extrapolate the linear evolution to the non-linear 
stage, but apply the Press-Schechter theory to interpret the resulting probability 
distribution. As for ^ 2 ^5 we estimate its growth rate by the linear theory, while 
we estimate its decay rate by including the effect of Compton drag. Specifically, 
we formulate the evolution of (^J 2 ^ as follows. 

First, we consider two models for the growth rate of the angular momentum. 



| A preliminary numerical evaluation of Pj suggests that its magnitude is enhanced at small 
|J| relative to the Gaussian case, which may affect our analysis below to certain extent. 
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Model 1 





(3.1; 



where 8 m = 1.06 which corresponds to the maximum expansion epoch of the spher- 
ical dust collapse (Fig. 2). Thus the linear theory is applied until this epoch in this 
model. 

Model 2 




f(j 2 )/t for 8 < 8 C , 



for 8 > 8 



(3.2) 



C 1 



where <5 C =1.69 which corresponds to the collapse epoch of the spherical dust collapse 
(Fig. 2). Hence, in this model, the linear growth rate of the angular momentum is 
extrapolated until the last moment of collapse if the cloud were a homogeneous 
dust sphere. 

On may regard the above two models as representing two extreme cases. In 
the former model, 8 = 8 m is roughly the validity limit of the linear theory. So, 
in this model, the growth of the angular momentum in the non-linear phase is 
ignored. Hence we expect it to give the final angular momentum less than the real 
value. At the non-linear stage, the growth rate will decrease due to the decrease 
in the quadrupole moment of the cloud. So, the latter model, in which the linear 
growth rate is extrapolated to the final stage of collapse, is expected to give the 
final angular momentum larger than the real value. 

Next, we consider the decay of the angular momentum due to Compton drag. 
We assume a sudden ionization history as 



. 10" 4 for z > 

Xe( z ) , n , 

1 lor z < z. 



where x e ^ s ^ ne ionization rate. With this ionization history, we model the decay 
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rate of the angular momentum as 




-2a Q (l + z) 4 {J 2 j for 8 < 1.69, 
for 8 > 1.69, 



(3.3) 



where 



a n\ z ) =%Xe{ z 




(3.4) 



7 n 



4cr T e 
3/xm c 





Here cr^ is the Thomson scattering cross section, e Q is the energy density of CMB 
today, m p is the proton mass, c is the speed of light, and // is the mean molecular 
weight which we assume to be 1, respectively. We cut the Compton drag at the 
collapse epoch, since the photon trapping occurs almost at this epoch (Loeb 1993). 

Now we solve the differential equation, 



(3.5) 



This can be analytically solved. We define the amplitude of 8 which characterize 
the reionization epoch, 

8. :=8 1 + Z . (3.6) 

ion 1+z. y J 

ion 

Then the solution is characterized by the value of 8. . We first show the results 

J ion 

for Model 1. 

Case 1A; 8. < 8 m . 

1 ion — 

In this case, the evolution of ( .T 2 \ is divided into the following four stages: 
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(i) 8 < 8. ■ 
P 



J2 

i 



l + z 



exp 



6a. t 



P((l + -) 5/2 -(l + ^) 5/2 ) 



J 2 



J 2 

ion 



1+2. 
wi 

l + Z 



exp 



6< Vo 



((i + 5/2 -(i + *) 5/2 )" 



(iii) 8 m < 8 < 8 C 



P 



■ J m 



exp 



((l + 5/2 -(l + -) 5/2 ) 



(iv) <5 C < 8 



P 



J 2 

■ J m 



const. 



(3.7) 



Case IB; <*L, < 8. < 8 n . 



We also have four different stages of the evolution in this case: 
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(i) 8<6 m ; 
P 



Jf) Kl + Z 



exp 



6a. t 



P((l + -) 5/2 -(l + .) 5/2 ) 



(ii) 8 m < 6 < 6. ; 

v / iii — ion 1 

6a. t 



P 



P 

■ J m 



exp 



p((l+^) 5/2 -(l+-) 5/2 ) 



(iii) 6. <6<6 C 

\ / inn — ^ 



P 



P 

ion 



exp 



6a 2 t Q 



((l + 5/2 -(l + -) 5/2 ) 



iv) S c < 6 ; 
P' 



P 

ion 



const.. 



Case 1C; 8„<8. . 



In this case, the evolution of ( P \ is divided into three stages: 



(i) S<6 m ; 



P 



P 



l + z 



exp 



6a. t 



P((l + ^.) 5/2 -(l + .) 5/2 ) 



(ii) 6 m < 6 < S c ; 



P 



J 2 

■ J m 



exp 



6a. t 



P((l + ^) 5/2 -(l + -) 5/2 ) 



(iii) S c < 6 ; 
P' 



P 

■ J m 



const.. 



In the above, the subscript "z" denotes the initial epoch which we take to be the 
decoupling time, "m" the epoch at which 8 = 8 m and "z'on" at z = z. . 

As for Model 2 } the evolutionary behavior becomes a bit simpler, since the 
epoch 8 = 8 m plays no role any more. 

Case 2A: 8. < 8 r . 

' ion — G 

The evolution of (<// is divided into three stages: 



(i) 8 < 8. ■ 



P 



J2 

i 



1 + z 



exp 



1 + zfl 2 -(1 + z 



,5/2 



(ii) 8. <8<8 C ; 



J< 



J2 

ion 



1+Z. 
101 

1 + Z 



exp 



6a 2 t Q 



1 + z ; 



,5/2 



:i + z 



,5/2 



(iii) 8 C < 8; 
P' 



const. 



(3.10) 



Case 2B; 6 m < 8. < 8 r . 



The evolution of ( J 2 ) is divided into two stages: 



(i)*<* c ; 



P 



J 2 



1 + z 



exp 



1 + zfl 2 -(1 + z 



>5/2 



(3.11) 



ii) 8 C < 8 ; 
P' 



const., 
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where the subscript "c" denotes the epoch 8 = 8 C . 

From the above results, the final value of the angular momentum after a clould 
is ready to collapse adiabatically can be estimated. For Model 1, we obtain 

Cases 1A and IB; 8. < 8 r : 

' ion — G 



J 1/ 2 -10/3 ( 1 + z i 
J2\ \l + Z f 



exp 



6t n 



^ a, (!+.-.)»/» + , )»/» 



(3.12) 



+ « 2 (( 1 + 2 .«) 5/2 -( 1 + - )5/2 )) 



Case 1C; 8 n < 8. 



J f 



-10/3 



J? 



1+^ 



(3.13) 



exp 



6a, t 



^ ((1 + -) 5/2 " (1 + - c ) 5/2 ) 



and for Model 2, we obtain 



Case 2A; 8. < 8 • 







! + ^- 


<*> 


■=( 




> 


s exp 


"_6^ 



^ a, (!+.-.)»/» + , )»/» 



(3.14) 



+ « S ((l + * j J ,/J -( 1 + --c>* /J )) 



Case 2B; 5, < 8. 
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where y^jj denotes the final mean square value of the angular momentum. We 
see that the final rms angular momentum for Model 2 is larger than that for Model 
1 by the factor 2 5 / 3 = 3.17. 

These solutions are plotted in Fig. 3 for a variety of z. . The amplitude of 
the density perturbation is set to the rms value estimated from the standard CDM 
model at scale M = 10 10 M Q . As easily seen from the figure, the final value of the 
angular momentum depends extremely sensitively on z. . As we have seen in the 
above, it also depends on the model of d ^ J 2 ^ / dt but not much compared to the 
z. -dependence. 

ion 1 

Having obtained the evolutionary model for ^ J* 2 ^ , we now have that of the joint 

— * 

probability distribution P (J , 8;t). Note that it does not factorize any more, due 
to the (^-dependence of (^J 2 ^ we have introduced. Using thus obtained probability 
distribution, we formulate a method to calculate the collasped mass fraction. As 
for 6, we use the Press-Schechter theory to estimate the collapsed fraction. As for 

— * — * 

J, we sum up the fractions which satisfy \J\ < J BH -> where J BH is the maximum 
angular momentum which a Kerr black hole can hold, 



GM h 

j bh = -t> ( 3 - 16 ) 



where M fe is the baryon mass contained in the mass scale M. 

Thus the mass fraction of baryons which collapsed into massive black holes on 
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scale M by the time t is 



d JP J8 (J,S;t 



\J\<Jbh 

oo 

dS 1 



; exp 



(3.17) 



Jbh 



J 2 dJ- 



47T 



2vr^J 2 



■ exp 



2^J 2 



4. Analysis and results 



In this section, we apply the formulae derived in the previous section to a 
cosmological model with the standard CDM density perturbation spectrum and 
present the results. We assume O fe = 0.05 and h = 0.5. 



Let us first consider the scale dependence of ( Jj 



the 



same as 



that of ( J? ) , which is evaluated 



as 



Since ( Jj) oc ( J() , it is 



■ J i 



1/2 



45 h 



p (ar 



,5 



0/ 



t 

15V4vry h (l + Zi ) 5 / 2 p 2 /H Q 



(4.i; 



In the above, the value of (8(t.) 2 ^ (or z c ) depends on the scale M in general, but 
within the range of our interest, 10 8 M Q ^ M & 10 10 M Q (corresponding to the 
black hole mass of 10 6 M Q £ M £ 10 8 M Q ), its dependence is very weak for the 
CDM spectrum. Hence we may conclude that 



j 



1/2 



OC 



(4.2) 



On the other hand, form Eq.(3.16), the maximum angular momentum of a black 
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hole is 

GM 2 

c 



j bh = . (4.3) 



Comparing this scale dependence of J Rff with that of ( J? ) , we find that there 



'BH VV1U11 W1 \ J f, 

exists a critical scale M cr at which ^Jf^j = ^bh anc ^ we nave ^ ^5_ff ^ or 

M > M cr . Hence density perturbations with the average amplitude can collapse 
to a black hole for M > M cr . This critical scale is calculated to be 



M 1 = _^_a-3^(i + ^)-3/2f% 

crl 32 • 125 6 c V c; G 



exp 



■^(«i((l + -.-) 5/2 -(l + 5/2 ) ( 4 - 4 ) 



^5/2 _ fU , W 2 



for Model L and M _ = 32M , for Model 2. They are plotted against z. in 

' cri crl ° ion 

Fig.4. 

The result for the black hole mass fraction is plotted in Fig. 5. As was 

the case of Fig. 3, the reionization epoch z. plays the most important role in 
determination of the final black hole mass fraction. For 700, more than 

ion 1 

10% of all the baryons will collapse into massive black holes. But for z. ^ 500, 

J 1 ion 1 

only an infinitesimal fraction will collapse into black holes. This result is almost 
independent of the model as well as of the mass scale. 



5. Summary 

We have calculated the joint probability distribution for the mean density con- 

— * 

trast 8 and the angular momentum J of a sphere from the Gaussian distribution of 
the density perturbation. We have found there is no correlation between the two 
variables in the linear regime. Then we have applied the Press-Schechter formalism 
to the density contrast and taken account of the Compton drag effect on the evo- 
lution of the angular momentum. By this procedure, we have introduced certain 

— * 

degrees of correlation betweeen 8 and J. 
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Using thus calculated probability distribution, we have evaluated the mass frac- 
tion of bound objects and distribution of their angular momentum. Then comparing 
this final angular momentum with the maximum possible one for a black hole, we 
have estimated the fraction which collapsed into black holes. We have found that 
the result is extremely sensitive to the epoch of reionization which we treated as 
a parameter. More quantitatively, compared with the observed number density of 
QSOs at redshift z ~ 2, the reionization epoch (z. ) must satisfy at least z. ^ 500 

^ ' 1 v ion' J ion 

(Fig. 5). We remark that this result is rather insensitive to the amplitude of density 
perturbations, because the acceptable range of reionization epoch is constrained by 
the Compton drag time scale. 

Of course, the present analysis is only a preliminary attempt to understand 
the history of the universe during the epoch 1000 & z & 10. For example, one 
future direction we should take is to perform a more self-consistent treatment of 
the reionization and the evolution of density perturbations. Also, it is necessary to 
calculate the probability distribution of the angular momentum in more detail and 
examine how much it differs from the Gaussian form, which we have assumed in 
this paper. 

ACKNOWLEDGEMENTS 

We would like to thank H. Sato, and R. Nishi for useful discussions. This work 
was supported in part by the Monbusho Grant-in- Aid for Scientific Research Fund, 
Nos.2010, 3077 and 05640342. 



- 19 - 



REFERENCES 

1. M. Kawasaki and M. Fukugita, preprint (1993). 

2. P.J.E. Peebles, Astrophys. J. 315 (1987), L73. 

3. K.M. Gorski, Astrophys. J. 398 (1992), L5. 

4. G. Efstathiou and M.J. Rees, Mon. Not. R. astro. Soc. 230 (1988), 5. 

5. M.G. Haehnelt and M.J. Rees, preprint (1993) 

6. W.H. Press and P. Schechter, Astrophys. J. 187 (1974), 425. 

7. Y. Suto, Prog. Theor. Phys. 90 (1993), 1173. 

8. P.J.E. Peebles, Astrophys. J. 155 (1969), 393. 

9. A.G. Droshkevich, Astrofisika 6 (1970), 581. 

10. S.D.M. White, Astrophys. J. 286 (1984), 38. 

11. B.S. Ryden, Astrophys. J. 329 (1988), 589. 

12. T. Quinn and J. Binney, Mon. Not. R. astro. Soc. 255 (1992), 729. 

13. A. Loeb, Astrophys. J. 403 (1993), 542. 

14. M. Umemura, A. Loeb and E.L. Turner, Submitted to Astrophys. J. Letters 
(1993). 

15. J.M. Bardeen, J.R. Bond, N. Kaiser and A.S. Szalay, Astrophys. J. 304 
(1986), 15. 

16. G.F. Smoot et al., Astrophys. J. 396 (1992), LI. 

17. J. Barnes and G. Efstathiou, Astrophys. J. 319 (1987), 575. 



- 20 - 



FIGURE CAPTIONS 



1. The existence of the wave vector k' corresponding to any wave vector k' , for 
a fixed k N , is shown. The contribution of the vector k' completely cancels 
that of k' in the sum of Eq.(2.17). 

2. A schematic picture of the radius of a collapsing homogeneous dust sphere as 
a function of the corresponding density contrast 8 in linear theory. The radius 
reaches its maximum at the epoch when the linear theory gives 8 = 1.06 and 
collapses to zero when 8 = 1.69. 

3. The evolution of the mean square of the angular momentum on scale M = 
10 10 M Q . Fig. 3a and Fig. 3b correspond to the two models in the text. In both 
figures, different curves correspond to different epochs of the reionization z. 
as indicated. 

4. The critical mass scale, below which a spherical overdensity region of of the 
rms perturbation amplitude can collapse to form a black hole, is plotted 
against the reionization epoch. 

5. The mass fraction of baryons which collapsed into massive black holes. The 
mass scale is M = 10 10 M Q , and we have adopted the CDM spectrum with 
h = 0.5 and = 1, normalized with the COBE DMR data at 10°. 

6. The same as Fig. 5 but for the mass scale M = 10 8 M . 
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